Nonlocality-induced front interaction enhancement 
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We demonstrate that nonlocal coupling strongly influences the dynamics of fronts connecting 
two equivalent states. In two prototype models we observe a large amplification in the interaction 
strength between two opposite fronts increasing front velocities several orders of magnitude. By 
analyzing the spatial dynamics we prove that way beyond quantitative effects, nonlocal terms can 
also change the overall qualitative picture by inducing oscillations in the front profile. This leads 
to a mechanism for the formation of localized structures not present for local interactions. Finally, 
nonlocal coupling can induce a steep broadening of localized structures, eventually annihilating 
them. 
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Most studies on the emergence of complex behavior 
in spatially extended systems consider that spatial cou- 
pling is either local or alternatively global (all to all cou- 
pling) [l|, More recently, systems with nonlocal (or 
intermediate- to long-range) coupling have received in- 
creasing attention, as nonlocal interactions are known to 
be relevant in diverse fields, ranging from Josephson junc- 
tion arrays ^ and chemical reactions 0, 0|, to several 
problems in Biology and Ecology @, such as the neural 
networks underlying moUusk patterns ocular dom- 

inance stripes and hallucination patterns 9|, and popu- 
lation dynamics [lol | . A nonlocal interaction may emerge 
from a physical/chemical mechanism that couples points 
far apart in space, e.g., a long-range interaction 11 11 , o r 
from the adiabatic elimination of a slow variable Il2l Il3| . 
Novel phenomena emerging genuinely from nonlocality, 
such as power-law correlations 12 j_ 3 , multiafHne turbu- 
lence Q , and chimera states |l5Lll6l| have been reported. 
Moreover, recent works have reported the effects of non- 
locality on the dynamics of fronts, patterns and localized 
structures (LS), for instance the tilting of snaking bifur- 
cation lines [17| and changes in the size of LS [ISj, the 
effects of two-point nonlocality on convective instabilities 
ligl, the nonlocal stabilization of vortex beams |20|, or 



changes in the interaction between solitons 21|, and in 



the velocity of propagating fronts [22 1 



The main goal of this Letter is to show the crucial 
relevance of nonlocality on the interaction of fronts con- 
necting two equivalent states in one dimensional systems, 
as well as on the formation of LS arising from the inter- 
action of two such fronts 23[. Interaction between two 



monotonic fronts is always attractive, so any domain of 
one state embedded in the other shrinks and disappears. 
However, fronts with oscillatory tails can lock at specific 
distances leading to stable LS. Here we show that oscil- 
latory tails, and therefore stable LS can appear as an 
effect of repulsive nonlocal interactions. Repulsive (in- 
hibitory) interactions are common, for instance, in neural 
field theories and genetic networks • Our result 



is generic and can be qualitatively understood from the 
interplay between nonlocality, which couples both sides 
of the front, and repulsiveness which induces a small de- 
pression at the lower side and a small hill at the upper 
part. Altogether this leads to oscillatory tails as qualita- 
tively obtained from the spatial dynamics. 

A prototypical model of a spatially extended system 
with two equivalent steady states is the real Ginzburg- 
Landau equation (GLE) [2|,|^ 
local real GLE 



25[- We consider the ID non- 



dE{x) 
dt 



(a* - s)E{x) + 



E%x)+sF{x), (1) 



being E{x) a real field and /i the gain parameter. The 
diffusion and the nonlinear term have been scaled to one 
without loss of generality. Eq. ([ij has both a local (dif- 
fusive) and a (linear) nonlocal spatial coupling 



F{x) 



r{x - x')E{x')dx' , 



(2) 



where 9^ is the spatial nonlocal interaction function (or 
kernel) and cr indicates the spatial extension of the cou- 
pling. For the sake of simplicity, we consider here a Gaus- 
sian kernel, O^ix - x') = (l/V2^cr)e-('"-="')'/(2'^'), al- 
though the results presented in this work do not depend 
qualitatively on its precise shape, provided it is positively 
defined. Gaussian kernels appear in contexts such as mol- 
lusk pigmentation patterns _26] and Neuroscience |§, . 
The nonlocal function F(x) includes also a local contri- 
bution. This contribution is compensated for by the term 
—sE{x), such that in the limit cr — one recovers the 
same results as for the local GLE. 

The profile of GLE fronts is known analytically p^. 
and the interaction with an opposite front located at a 
distance d much larger than the front width can be cal- 
culated perturbatively, yielding the following result for 
the relative velocity v{d) [28] . 



v{d) = d — ce 
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FIG. 1: Front velocity as a function of d. ^ = 3, |s| = 1, and 
a is taken to be , 2 and 10, depicted in the figure as cr''^"'^'. 




FIG. 3: (Color online) Nonlocal GLE, Eq. ^ with /i = 3: 
(a) Front profiles for s = Q (black line) and s = ±1 with 
o = 2, labeled in the figure as cr=s"'"'. (b) The curve A 
is the boundary in the (cr, s) space between the presence of 
oscillatory or monotonic tails of a front. 
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FIG. 2: 7 as a function of cr. Inset: 7 as a function of 1 /a. The 
black solid, dark gray dashed and light gray dashed curves 
have = 3, and s = 0.5, 1, and 2, respectively. 



with c = —24-^/2/7 and 7 = V2/i. The nonlocal effects in 
front dynamics can be quantified by looking at the devi- 
ations from the interaction given by Eq. ([3]), stemming 
solely from a local interaction coupling. Therefore we 
have studied the front velocity v{d) for different kernel 
widths fj keeping all other parameters fixed and taking 
the system size much larger than cr. The results are plot- 
ted in Fig. [TJ The local case, Eq. ([S]), is the straight line 
labeled with (as cr = 0). In the nonlocal case (cr 7^ 0), 
a large qualitative difference in behavior is observed de- 
pending on the sign of s. 

As shown in Fig. [1] for an attractive (activatory) in- 
teraction (s > 0), the logarithm of the velocity decreases 
linearly with the distance, such that the exponential de- 
pendence of the velocity with the distance given by Eq. 
^ still holds with an effective 7 whose value is strongly 
reduced by the nonlocality. As a consequence the range of 
spatial interaction, I/7, increases with the kernel width 
cr. Even for moderate values of ct, fronts move several 
orders of magnitude faster than for the local GLE. Fig. 
[5] shows the change of the value of the effective exponent 
7 as a function of cr. Moreover, the inset shows that 7, 



to a very good extent, follows a linear dependence with 
the inverse of the kernel width cr, provided that cr is at 
least as large as the front width. Therefore, we can con- 
clude that rescaling 7 to 7/CT the effective interaction 
between two fronts follows a universal exponential law 
(except for a small dependence on the strength of inter- 
action s). It is interesting to notice, however, that the 
width of the front, defined by the half width at half the 
maximum (HWHM) does not show this scaling with l/cr 
cf. Fig. ^a). This can be understood by noticing that 
while the general shape of the front is mainly dominated 
by the local, diffusive, coupling, in turn the nonlocal cou- 
pling modifies substantially the exponential tails, which 
are responsible of the long-range interaction that influ- 
ences the front velocity. 

For s < 0, i.e., for a repulsive (inhibitory) interaction, 
the exponential law, Eq. ([3]), no longer holds. Neverthe- 
less, the magnitude of the envelope of the front velocity 
is still dominated qualitatively by ([3]), as shown in Fig. 
[T] In this case, the velocity becomes zero at regular in- 
tervals of the distance d between two fronts. At these 
positions the fronts are locked leading to the formation 
of LS. These LS, not present in the GLE, come into exis- 
tence by the creation of oscillatory tails when a nonlocal 
term with s < is added. 

A quantitative characterization is obtained from the 
spatial dynamics of Eq. ([T]). Considering a spatial per- 
turbation to the homogeneous solution of the form E = 
^/JI + eexp(Aa;), one finds 

-2^i-s + X'^ + se^''"'/2 ^ 0. (4) 

Eigenvalues come in pairs ±A. The shape of the front 
is determined by the eigenvalue Ai with real part clos- 
est to zero, as in the spatial dynamics all the other di- 
rections are damped faster when approaching the fixed 
point (homogeneous solution). By determining where Ai 
goes from being purely real to complex, one can find the 
boundary A separating fronts having monotonic and os- 
cillatory tails. As shown in Fig. EJb), for kernel widths 
cr at least as large as the front width, the front profile 
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FIG. 4: Stationary widths of stable LS for the nonlocal GLE 
(/i = 3). Black (gray) line corresponds to a = 0.5 (cr = 
2). Insets show examples of LS profiles for (7 = 2. Scross,<y, 
given by line A in Fig. |31 indicates the point where the spatial 
eigenvalues go from being complex conjugate to purely real. 
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FIG. 5: Stationary widths of stable LS for a — 1, = 0, 
fj, = Q, V = 2 and p = 2.7. a = 1 (black line) and 2 (gray 
line). Insets show the transverse profile of the LS. 



always shows oscillatory tails for s < and as a conse- 
quence LS can arise as displayed in Fig. |4l For smaller 
kernel widths an increasingly large nonlocal strength is 
needed for LS to be formed (see line A in Fig. [3]). To 
describe the interaction of fronts with oscillatory tails 
Eq. ([2]) must be modified. An appropriate ansatz is the 
following: 



V — d ~ ccos.[C,{s,a)d\e 



(5) 



where C(s,a) is determined by the complex part of the 
most underdamped spatial eigenvalue of Eq. (|3]), such 
that C = in the region of monotonic tails. Eq. (O 
adequately describes the dependence of the velocity with 
the distance indicated in Fig. [TJ 

In summary, for the real GLE one finds an enhance- 
ment of the interaction due to nonlocality, that leads to 
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FIG. 6; (Color online) (a) Dependence of the absolute value 
of the real and imaginary parts of the two dominant complex 
quartets of spatial eigenvalues in the nonlocal PCGLE for 
(7 = 1.3 and cr = 2 as a function of s. Other parameters as in 
Fig. [5] (b) The curve B depicts the locus in the (a, s) plane 
of the points where the absolute values of the real parts of A; 
and \ni cross. 



the appearance of oscillatory tails and LS for a repul- 
sive interaction and increase of the front velocity if the 
interaction is attractive. 

We check the generality of these findings for the Para- 
metrically forced Complex Ginzburg Landau Equation 
(PCGLE) whose fronts show oscillatory tails already in 



the local case and form LS 



The PCGLE is the 



generic amplitude equation for oscillatory systems para- 
metrically forced at twice the natural frequency |31|] . The 
nonlocal version of the PCGLE can be written as: 



dE{x) 
dt 



= (1 



']E{x) 



~{l+iP)\E{x)\'' E{x)+pE*{x)+se'^'F{x), (6) 

where /i measures the distance from the oscillatory in- 
stability threshold, v is the detuning between the driving 
and the natural frequencies, p > is the forcing ampli- 
tude, a and (3 represent the linear and nonlinear disper- 
sion. The term se^^'F{x) describes the nonlocal response 
of the material taking the same form as in the study of 
the GLE (0). Here we consider s > with 0s = and 
(ps — TT in correspondence with attractive and repulsive 
interactions discussed previously. Again, the linear con- 
tribution has been compensated in the term —se^'^''E{x). 

Fig. [5] shows the width of the stable LS locked at the 
first and second oscillations as a function of the nonlo- 
cal strength s. The insets show the LS spatial profile. 
Whereas in the real GLE, LS existed only for s < (see 
Fig. U), in the PCGLE, LS also exist for a finite range 
of positive values of scos((/)s) {4>s = 0). This is shown in 
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Fig. [S] for (T = 2 where the bifurcation branch abruptly 
ends around Scross = 0.556. This can be explained again 
by studying the spatial eigenvalues of the system. In this 
case there are two pairs of complex conjugate eigenval- 
ues with small real part A; and A„; which play a relevant 
role. Fig. IHl^a) shows the dependence of the real and 
imaginary parts of these two eigenvalues for a — 1.3 and 
2 as a function of s. For small s the spatial dynamics 
is dominated by A;, an eigenvalue already present in the 
local case. Increasing s the real part of this eigenvalue 
becomes slightly larger while the real part of Xni clearly 
decreases. For a — 2, the real parts of the two eigenval- 
ues cross at Scross = 0.556. Beyond that value the spatial 
dynamics is governed by A„; whose spatial length scale 
[Im(A„;)] is an order of magnitude larger than the one 
corresponding to A;. As a result, the branch of LS expe- 
riences dramatic sharpening in Fig. [SJ as the LS broaden 
with an order of magnitude. Fig. Hljb) shows the locus 
of the crossing point in the (cr, s) plane for ips — 0. For 
large values of a, the crossing point moves towards s — 0, 
while for a < 1.43 the crossing never takes place and the 
spatial dynamics is always dominated by A;. 

For = IT., nonlocality only modifies slightly the pre- 
existing oscillatory tails of the PCGLE. The branches of 
LS end for negative values of s cos(0s) close to the mod- 
ulational instability point of the background states (that 
for cr = 2 is at scos(0s) ^ —3). This MI point, that 
was originated by the nonlocality in the GLE, is already 
present in the local form of the PCGLE but is enhanced 
by the nonlocal interactions. 

In conclusion, we have demonstrated the large impact 
of a linear nonlocal term on the interaction of fronts con- 
necting two equivalent states. The most striking result is 
the possibility of obtaining a novel class of self-organized 
stable localized structures in systems exhibiting fronts 
with no tails, like the real GLE. This is achieved through 
a nonlocal repulsive interaction that modifies the profile 
of the fronts introducing or damping out oscillations. In 
addition, we have observed an order of magnitude in- 
crease of the front velocity due to an enhancement of the 
interaction between two fronts. Finally we have shown 
that nonlocal interactions can also smooth out the front 
oscillations greatly increasing their wavelength such that 
the LS become much wider and eventually disappear. 

The characterization of the effects of nonlocal coupling 
on the front properties and dynamics can allow the iden- 
tification, from both theoretical and experimental data, 
of different sources of nonlocality. For instance, domain 
walls have long been studied in photorefractive media 
(3^ , which have a large nonlocal response, but its effects 
on the front dynamics have never been identified. Nonlo- 
cal interactions are also common in Biology, Chemistry, 
and other fields, and they can have a constructive role by 
enhancing the propagation of information between dis- 
tant parts of the system. 
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